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Abstract 

We point out that the matter profile effect — the effect of matter density fluctu- 
ation on the baseline — is very important to estimate the parameters in a neutrino 
factory with a very long baseline. To make it clear, we propose the method of the 
Fourier series expansion of the matter profile. By using this method, we can take 
account of both the matter profile effect and its ambiguity. For very long baseline 
experiment, such as L = 7332km, in the analysis of the oscillation phenomena we 
need to introduce a new parameter ai — the Fourier coefficient of the matter profile 
— as a theoretical parameter to deal with the matter profile effects. 



1 Introduction 



The evidence of the neutrino oscillation was presented by the Super-Kamiokande collaboration [|I|]. 
Many experiments have been planed and some of them have already been carried out to 
determine the parameters of the neutrino oscillation0. Recently, the neutrino factory]^ 
— a very long (over 1000 km) baseline experiment with high energy neutrinos from a 
muon storage ring — has been proposed as the most effective experiment to determine 
the unknown parameters, Ues, the sign of Am^ which is responsible for the atmospheric 
neutrino anomaly and the CP violating phase |, |, |, 0, g ^ 0, |Tl|, 0, |T3[ . 

To analyze the oscillation phenomena with such a long baseline, we have to treat the 
matter effect very carefully. Two ways have been adopted to deal with the matter 
effects: (i) To add the averaged matter density p as a parameter in addition to the other 
theoretical parameters like the mixing angles, and estimate it in the same way as the other 
parameters[^, ^, |ll|]. (ii) To use the PREM — Preliminary Reference Earth Model[15 — 



for the matter density profile, and to assume that we know the ingredients of the earth 
completely!^, ^f. We have, however, questions to these treatments. For the method-(i), 
"Can we describe the matter effect precisely enough by only averaged matter density 
?" The answer was "No, we can't. Sometimes the deviation from the constant density is 
important. " 0]. For the method- (ii), "Is the PREM a trustworthy model for the neutrino 
oscillation experiments ?" The PREM is originally based on the study of earthquake 
waves, hence it predicts the density profile in the depth However, it doesn't predict the 
ingredients of the matter. Therefore, there is an ambiguity in electron number fraction. 
Then, we should worry whether there are ambiguities in electron number density. If we 
estimated the parameters without considering its ambiguity, we would have significant 
errors in the estimates for the parameters. Therefore, we must introduce the way to take 
into account the matter profile effect as the parameter determined by the experiments. 
We propose the method of the Fourier series expansion of the matter profile p6[]. Using 
the Fourier expansion we can parameterize the matter profile. We can express the matter 
profile effect with a finite number of the parameters by examining how many terms of 
the Fourier expansion contribute to the oscillation physics within the resolution of the 
experiments. We can incorporate the matter profile's ambiguity in the ambiguities of the 
Fourier coefficients. 

We will review the method of the Fourier series and investigate qualitatively what 
circumstances make the matter profile effect relevant in section 2. In section 3, we will 
calculate quantitatively the oscillation probabilities and the event rates with various sets 
of the parameters in baseline lengths. 



2 The method of the Fourier series 

In this section we introduce the Fourier expansion method. First, we describe its formal- 
ism and see its feature. Then, we solve the evolution equation for neutrinos in matter 
perturbatively, and study the condition where the matter profile effect is significant. 
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2.1 Fourier expansion method 

Assuming three generations of neutrinos, we parameterize the lepton mixing matrix 

'-^ m — ^ <-^Majorana 




^Majorana 

1/ 

CtpSi^ S^S^Ci^C C-ipC^jj SipS(j)S,jj6 SipC(j)C I f^Majorana- (1) 

which relates the flavor eigenstates = e, /x, r) with the mass eigenstates in vacuum 

= 1, 2, 3) as 

Wa) = UaiWi). (2) 

^Majorana IS the part of the Majorana phases. It doesn't contribute to the neutrino oscil- 
lation phenomena and hence is omitted hereafter. 

The evolution in matter of the flavor eigenstates of neutrinosQ with its energy E is 
given by 



. d 

dx 






fa{x) ^ 














f3aJ 





where 



lypix)) = H{x)f3aMx)) (3) 



Am^^- = rrij — , mi : mass eigenvalue 
a(.) . 2V-2G.n.(.)E = 7.56 x 10- (^) (^) K, (4) 



ne{x) 
p{x) 



Fermi constant 
electron number density 
matter density 
electron fraction 



To parameterize the matter effect, a(x), we expand a{x) into Fourier series as 

°° 27r 
«(^) = XI = -j^n. (5) 



^For antincutrino a{x) and S should be replaced by —a{x) and —S. 
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Note that a_„ = a* due to the reahty of a{x). Also within the PREM, 

a-n = an (6) 

since a{x) = a{L — x). 

If this expansion can be approximated with a finite number (= A^) of terms, 

n=N 

a{x) = «ne-*^"^ (7) 

n=-N 

then it means that the matter profile effect can be parametrized. Thus, by introducing 
new parameters a„(n = — A^, ■ ■ ■ , A^) (orn = 0, ■ ■ ■ , A^, if eq.(|D holds), we can investigate 
the oscillation physics without help of earth models. 



2.2 Perturbative analysis of matter profile eff'ect 

We solve the evolution equation (0) to see the qualitative feature of the matter profile 
effect. First, we divide the Hamiltonian H{x) into three pieces for the later calculation: 



00 







— e'^^'T Am2icl I T^e-'^^\ 

^Aml^c^s^ - Amj^c^s^sl Amj^cl + Amais^s^ 



-f^oi = 



2E 







(9) 



Hiix) = — 



6a(x) 



0, 



(10) 



where a = Oq and 6a{x) = a{x) — a. The oscillation phenomena arise from Hqq, hence its 
effect must be fully taken into account. On the other hand Hqi and Hi{x) can be treated 
as the perturbation in almost all cases in neutrino factories, since 



Aml^L 



2E 



< 1 and 



6a{x)L 



2E 



< 1 



can be almost always expected. Note that we have to include Am^i terms in {HQQ)ij{i,j 
1,3) to deal with the resonance effect .0 
Next, we rewrite Hqq as 



1 ~ z^'- 

^^00 = ^Uo I Aml.cl 







(12) 



^On the other hand, we include Am^i term in the (22) element just for convenience. This inclusion 
is not essential. 
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where Lfo and A± are the mixing matrix and mass square eigenvalues in matter up to the 
zeroth order of the "perturbation" s Hqi and Hi{x) and defined by 



C2 







tan 20 



(13) 

(14) 
(15) 



a 



2 2 
21 ■'a; 



/3 ^ J {{Ami, - Amlsl)c2^ - ay + {A 



Then we solve the evolution equation (|^) perturbatively as 

-L 



exp 



H{x)dx 



Km 



I3a 



Hoijdx + e 

(5*00 + 5*01 + Si H )/3a|z^a(0)), 



-iHooL I 



Hi{x)jdx + 



Pa 



Km 

(16) 



where Sqo, •S'oi and 5*1 are defined by the first, second and third terms in the second last 
line respectively. 5*01 and 5*1 are the perturbative contributions at the first order from Hqi 
and Hi respectively. 5*00 , •S'oi and Si are calculated to be 



'S'oo(-^)/3a 



-iHooL 



Un 



( 



\ 



. Am 



2 2 
21^L 



1A5 



e -2E^ — 



I AS 



e 2B — e 2B 



|s2,/,(c?^e"*2B'^ + s|e~*2B^ — e" 



1 ^-j<5 



e 2B — e 2B ' 



^■52^/)(c^e"*2B'^ + s|e~*2B-^ — e"* 2b -^) 
c^(c^e"*2E'^ + s4e"*2B^) + s^e~* 2b y 



(17) 



'S'oi(-^^)/3a = e 



-I] 



dxHi 



01/ 



A 0^ 
Uo\A i3 I f/J 
B 



-s^{cgA + sgB) 



c^{c^A + s^B) 
-S2i,{s^A-c^B)cs 
[e^^cl-e-^'sl){s;i,A-c^B) 



-s^{ci,A+ s^B) 
{e^'sl-e-%)is^A-c^B) 
S2^is^A- c^B)cs 

(18) 



A 



B = 



2{Amlcl - A_ 



^ OZ7' ^''OIT' ^ 



e " 2B 



2(A+ - Aml.cl 



(19) 
(20) 



Cn 

E^o ( I t/o^ 




^ e^'5s^(P„c^ + C„sp -(C„ + r'„)4s^c^ -(C„ + P„)s^cv,s^c^ 
n \^e%(P„c^ + C„sp -(C„ + r'„)s^c^s^c^ -(C„ + r'„)4s^c^ y 



E 



S20(e~*2E — e~*2B''") 



(21) 



"2{(A+-A_)2-(2i5;p„)2} 

(A+ - A_)s2^ e-^-^s^ {(A+ - A_)c2^ - 2Ep„} 6"% {(A+ - A_)c2^ - 2Ep4^ 

X I e'^s^, {(A+ - A_)c2^ + 2Epn} -(A+ - A_)4s2^ |(A+ - A_)s2v,s. 

^e%{(A+-A_)c2^ + 2Ep4 



|(A+ - \-)S2^i 



2<f) 



-(A+ - A_)C^S20 



C = 



2(A+ - A„ + 2Epn] 

On 

2(A+ - A_ - 2Epn 



'^24> I ^ — e '^zB- 



(22) 
(23) 



With eq.(|l7|), (|18]) and (^T]) , the oscillation probability, for example from to i/^, up 
to the lowest order of the perturbations, Hqi and Hi{x), is calculated as 



P.^^. = |^(L) 



fie I 



|(5oo);.er + 2i?e[(S'oo);.e(5oi):e ] + 2i?e [(5'oo);.e (^1 ) 



/ fiei 



5 



2 I 

AE 



+ 



ClCj 



A"^21 



+ ^S5S2V.S2a;S2^ 



- A_ 
Am^^ 



SlS 







- Aml^clJ 








- - A"^2icS 




Am^^ 



Am^,c5-A,^ 
- Am2,c2 



sm 



'Am2^c2 



A 



Am: 



A. 



- A_ 



L 



sm 



'21 "-a; 

Am^^ 



4E 



X sm 



A, 



Am^i cl ^ ^ . Am^i c; 
-L + sm - 



2 2 
2lCij 



A 



Am2^c2 



(A+ - A_)C2^ + 2Epn 



2E 



A. 



L- 
A_ 



A+-A, 



(Ah 



A_l 



i2Epn) 



■ sm 



sm 



-L. 



2E 



(24) 



From this equation we find the following; (i) The matter profile effect is relevant when 
some of the Fourier coefficients large as Am^i or the resonance condition |TT 



A+-A_ 



2Ep„ 



(25) 



is satisfied, (ii) The matter profile effect decreases in proportion to 1/n when the resonance 
condition is not satisfied. Therefore the higher Fourier modes are expected to be irrelevant. 
In other words, we don't need the detail of the matter profile. Thus we expect that we 
can truncate the Fourier expansion of a{x) at the finite number N. We can estimate the 
matter profile effect in neutrino factories. 



3 Numerical analysis 

In this section we examine how many terms in the Fourier expansion are necessary. To see 
this we compare the oscillation probability and the event rate calculated using the PREM 
with those calculated using the truncated Fourier series (|^. We first calculate the transi- 
tion probabilities with several matter profiles, the PREM matter profile and the Fourier 
series matter profiles that are truncated at various iV's. Next using these probabilities 
we derive the event rates in the appearance channel, — > z/^ as the "observed" numbers 
in an experiment. To derive them we use the following equation with a given number of 
muons in decay A^^, muon energy E^, detector size N^t and detection efficiency e. 

flux 

' •> probability detector 
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ay = 0.67 X 10-^^ X (^-^y^ [m^] : CC cross section 

A^^ : Avogadro's number 
_ 

7 = , m„ : muon mass 

Eth '■ threshold energy of the detection. 

Then we compare these results. Through these comparisons we find how many Fourier 
coefficients we have to introduce to the theoretical parameters. 

In the following we assume that x Ny^t x e equals to = 5 x 10^^ kt.0 

3.1 L = 3000km 

Neutrino factories of this length are considered to be most efficient, so we consider this 
case first. 

The matter profile^ and its Fourier coefficients are shown in Fig.|I[ From these figures 
we find that the density fiuctuation is very small. It should also be noted that the 
resonance condition ( ]25|) will not be satisfied for any n in neutrino factories. These facts 
indicate that the the constant matter density approximation is valid. In other words, the 
Fourier modes will be irrelevant. To confirm this we plot the oscillation probabilities in 
Fig.|| with the various matter profiles; (i) constant density (dotted line), (ii) eq.(|^) with 
N=l (dashed line), and (iii) PREM matter profile in Fig.|l| (solid line). We find that 
indeed this figure confirms our expectation. 

The event rates are shown in Fig.|^. We plot the event rates calculated with the 
following 5's and the matter profiles; (i) 6 = 7i/2 and constant density (dotted fine), (ii) 
S = 7r/2 and the PREM matter profile (solid line), and (iii) 5 = and the PREM matter 
profile (dash-dotted line). The event rate with 5 = 7r/2 is different from that with 5 = 0. 
We can observe this difference, if we know all the other theoretical parameters, the mixing 
angles, Am^'s and d accurately enough]^. Therefore we have to consider the contribution 
of S to the oscillation probability. However, we cannot see the difference of the event rate 
in the case that oi = and that ai ^ with baseline L = 3000km. Thus we do not have 
to take into account the ffuctuation of the matter density. 



3.2 L = 7332km 

There are many analyses of the oscillation physics with the baseline length L = 7332km 
too. Therefore we consider this case next. 

■^We also assume that we can reconstruct neutrino energy. 
To calculate the matter profile, we linearly interpolate the PREM's data given discretely for the 
depth. For this calculation we modify the density of the crust(= most outside layer) into 2.3 g/cm"^ from 
1 g/cm^. 

^We cannot see this difference, however, since we will not determine the theoretical parameters accu- 
rately enough 
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Figure 1: (a) Matter profile on the baseline and (b) its Fourier coefficients for L — 3000km. 
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[GeV] 

Figure 2: Transition probability Pu^^Uf, for L = 3000km. In this pfot we set sinc<j = 
l/V2,sinV^ = l/\/2,sin0 = 0.1, Am^^ = 3 x lO'^eV^ Am^^ = 5 x lO'^eV^ and 6 = 7i/2. 
There are three hnes corresponding to the matter profiles; (i) constant density (dotted 
line), (ii) eq.(|^) with N=l (dashed line), and (iii) PREM matter profile in Fig.|I| (solid 
line). These lines are quite similar to each other so that we cannot see those three lines 
in the graph. 
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Events per GeV 




(a) = 30 GeV 



[GeV] 




[GeV] 



Figure 3: Differential event rate N^^^i,^ for L = 3000km with (a) Efj_ = 30 GeV and (b) 
Efj_ = 50 GeV. The oscillation parameters are same as those in Fig.|] except 6. There 
are three lines with respect to 6 and the matter profiles; (i) 5 = n/2 and constant density 
(dotted line), (ii) 5 = 7r/2 and PREM matter profile (solid line), and (iii) 5 = and 
PREM matter profile (dash-dotted line). The line (iii) is different from the others. 
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The matter profile based on the PREM and its Fourier coefficients are shown in Fig.|^. 
We find that the first two coefficients are large. We argue, however, that the first Fourier 
coefficient among them contributes much more to the oscillation probability because the 
resonance condition (|25|) for n = ±1 can be satisfied in a neutrino factory. 

In Fig.^, we compare the oscillation probabilities calculated with different matter 
profiles; (i) constant density (dotted line), (ii) eq.(^ with = 1 (dashed line), (iii) eq.(|^ 
with N = 2 (dash-dotted line), and (iv) the PREM matter profile in Fig.§ (solid line). 

We see that the probability calculated using PREM (Fig.^ - (iv)) and another cal- 
culated using eq.(0) with = 1 (Fig.^ - (ii)) are very similar. However the probability 
calculated using the constant density matter profile (Fig.^- (i)) is quite different. Fur- 
thermore there is no significant difference between the probabilities calculated using eq. (|^ 
with = 1 and N = 2 (Fig.|^ - (ii) and (iii)). Thus we have to take into account a new 
parameter ai in the analysis of neutrino factories with the baseline L = 7332km. This 
new parameter, ai, should be estimated by the experimental results. 

To see the importance of ai more clearly, we show the event rate in Figs.|^ using 
different 6 and different matter profiles; (i) 5 = 7r/2 and constant density approximation 
(dotted line), (ii) 5 = 7r/2 and the matter profile (0) with = 1 (dashed line), (iii) 
5 = tt/2 and the PREM matter profile (solid line), and (iv) 5 = and the PREM matter 
profile (dash-dotted line). 

The distributions of the event rate in (Figs.P-(i)) and (Figs.|]-(iii)) are quite different 
from each other. This means that the estimation of the event rate with constant density 
approximation doesn't work. In other words we cannot estimate the values of the param- 
eters from the experimental data precisely without taking into account the contribution 
of ai. On the other hand, the event rates in (ii) and (iii) are quite similar. If we can 
insist from Figs.^that we can observe CP violating effect, then we can insist from Figs.^ 
that we can observe the effect of ai. Furthermore since the contribution of ai cannot be 
explained by other terms, we expect that we can estimate ai very well experimentally. 
We can study geophysics by neutrinos. 

Consider the difference between Fig.^(iii) and Fig.p-(iv). The effect of the CP vio- 
lating phase on the event rate is so small that we cannot distinguish between the two 
theories that have different CP violating phase. Indeed the effect of 5 is smaller than that 
of 02 which is already beyond the experimental sensitivity. 



3.3 L = 12000km 

In this case neutrinos penetrate the earth almost along its diameter. Neutrinos go through 
both the mantle and the core. Therefore the matter profile is very complicated. We expect 
that we can observe the higher modes in eq.(0). 

The matter profile and the Fourier coefficients are given in Fig.|^. From these figures 
we expect that not only ai but also the higher Fourier modes are relevant. 

In Fig.|^, we plot the transition probabilities corresponding to the matter profiles; (i) 
constant density (dotted line), (ii) eq.(|^) with N=l (dashed line), (iii) eq.(^ with N=3 
(dash-dotted line), and (iv) PREM matter profile in Fig.^ (solid line). From FigJ^ we 
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Figure 4: (a) Matter profile on the baseline and (b) its Fourier coefficients for L = 7332km. 

12 




[GeV] 

Figure 5: Transition probability P^^^u^ for L = 7332km. The oscillation parameters are 
the same as those in Fig.^. There are four lines corresponding to the matter profiles; (i) 
constant density (dotted line), (ii) eq.(|^) with N=l (dashed line), (iii) eq.(^ with N=2 
(dash-dotted line), and (iv) PREM matter profile in Fig.^ (solid line). The line (i) is quite 
different form the others while the other three lines show quite similar line shape. 
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Figure 6: Differential event rate N^^^^,^ for L = 7332km witli (a) E"^ = 30 GeV and (b) 
Efj_ = 50 GeV. Tlie oscillation parameters are the same as those in Fig.^ except 6. There 
are four lines with respect to 6 and the matter profiles; (i) 6 = 7t/2 and constant density 
( dotted line), (ii) 5 = tt/2 and eq.(0) with N=l (dashed line), (iii) 6 = tt/2 and PREM 
matter profile (solid line), and (iv) 6 = and PREM matter profile (dash-dotted line). 
The line (i) is different from the others while the other three are quite similar. 
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Figure 7: (a) Matter profile on the baseline and (b) its Fourier coefficients for L — 
12000km. 1 




[GeV] 

Figure 8: Transition probability Pu^-,v^ for L = 12000km. In this plot we use the same 
oscillation parameters as those in Fig.^. There are four lines corresponding to the matter 
profiles; (i) constant density (dotted line), (ii) eq.(|^) with N=l (dashed line), (iii) eq.(|^ 
with N=3 (dash-dotted line), and (iv) PREM matter profile in Fig.0 (solid line). The 
line (i) is quite different from the others. The line (ii) is also different from the line (iV), 
while the line (iii) shows a quite similar shape with the line (iv). 
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Figure 9: Differential event rate N^^_^i,^ for L = 12000km with (a) = 30 GeV and (b) 
Efj^ = 50 GeV. The oscillation parameters are the same as those in Fig.^ except 6. There 
are three lines corresponding to 5's and the matter profiles; (i) S = 7r/2 and eq.(|^ with 
= 1 (dotted line), (ii) 6 = n/2 and eq.(0) with = 3 (dashed line), (iii) 6 = 7i/2 and 
PREM matter profile (solid line); (iv) 6 = and PREM matter profile (dash-dotted line). 
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find that there are significant contributions to the oscillation probability from the higher 
modes. The probability calculated with the constant density matter profile (Fig.H-(i)) 
is quite different from that calculated with the PREM matter profile (Fig.|^(iv)). Even 
the probability calculated using the matter profile eq.(0) with N=l (Fig.|^(ii)) differs 
apparently from that with the PREM matter profile (Fig.^(iv)). On the other hand we 
find that the matter profile eq.(|^) with N=3 mimic very well the PREM matter profile to 
calculate the oscillation probability. 

We also plot the event rates in Figs.^. There the lines correspond to the following 5 
and matter profiles; (i) 5 = n/2 eq.(0) with = 1 (dotted line), (ii) 6 = tt/2 and eq.(0) 
with = 3 (dashed line), (iii) 6 = tt/2 and PREM matter profile (solid line); (iv) 6 = 
and PREM matter profile (dash-dotted hne). We see little discrepancy in Figs.p-(ii), (iii), 
and (iv), while there is a conceivable difference between Figs.^(i) and the others. We 
find, therefore that the effect of 6 is irrelevant with the oscillation physics while the first 
three modes of the Fourier coefficients have measurable contribution to it. 

We examine more carefully whether the higher mode contribution is really measurable. 
The number of the necessary new parameters is quite dependent on the event rate of the 
experiment. We can find that there is conceivable difference between the lines of Figs.p-(i) 
and those of Figs.H-(ii)& (iii). However this difference is not significant statistically in 
Fig.^] (a). As long as A"^ x A^kt x e is less than 5 x 10^^, we need to introduce only ai as 
the theoretical parameter. On the other hand this discrepancy is very significant in Fig.^ 
(b). We can measure the contributions from the higher modes such as a2 and as. We 
should introduce higher modes according to the expected event rate. 

4 Summary and discussion 

We considered how to deal with the matter profile effect. We proposed that the Fourier 
coefficients of the matter profile are used as the theoretical parameters. Using this method 
we can evaluate the size of the ambiguities in the estimate of the mixing parameters. 

The perturbative solution for the evolution equation, eq.(^), implies that the higher 
Fourier modes are irrelevant. The introduction of the first few modes gives enough preci- 
sion to the estimate of the event rate. 

We saw the following three cases in detail numerically. 

• L = 3000km. The matter profile effect itself is irrelevant. We can assume the matter 
density is constant. 

• L = 7332km. We need to introduce ai as the theoretical parameter which should be 
measured experimentally. On the other hand the CP violating phase 6 is irrelevant. 
We cannot determine the CP violating phase at this baseline. 

• L = 12000km. We need to introduce higher modes as the theoretical parameters. 
However, since the event rate with this baseline length is significantly small, by using 
only oi we can estimate the theoretical parameters accurately within the precision 
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of the experiment. Depending on the precision required, we must decide the number 
of coefficients to introduce. 

In general the Fourier modes are relevant ,when the baseline passes through the lower 
mantle(L > 5000km). In this case we need to use the method developed in this paper 
to make an accurate analysis. In other words we can explore the interior of the earth in 
neutrino factories. 
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